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Masses for the (Qs)^*' (Qsy*' (Q = c or b) molecular states are systematically computed in the 
framework of QCD sum rules. Technically, contributions of the operators up to dimension six are 
included in operator product expansion (OPE). The numerical result 4.13 ± 0.10 GeV for D*D* 
agrees well with the mass 4143.0 ± 2.9 ± 1.2 MeV for F(4140), which supports the D*D* molecular 
configuration for F(4140). 

"q^" 1 PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 

O ' 

o 

I. INTRODUCTION 

Recently, the CDF Collaboration has reported the observation of a narrow near-threshold structure in 
the mass spectrum in B + — > J/%jxj)K + decays [l|, for which the mass is 4143.0 ± 2.9 ± 1.2 MeV 

and the width is II.T^qq ± 3.7 MeV. This experimental observation has triggered great interest of many 
' practitioners, and there have already appeared some theoretical interpretations for this new resonance, 
^ e.g. Refs. 0, E0, S, fl. On the whole, Y(4140) is apt to be deciphered as the molecular partner of 
the charmoniumdike state V(3930) 0]- Undoubtedly the quantitative description of V(4140)'s properties 
such as mass is quite needed for well understanding its structure, but it is difficult to extract information 
on the hadronic spectrum from the rather simple Lagrangian of QCD. That's because low energy QCD 
£SJ , involves a regime where it is futile to attempt perturbative calculations and one has to treat a genuinely 
str ong field in nonperturbative methods. Whereas, one can resort to QCD sum rules Q (for reviews see 
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12| and references therein), which is a nonperturbative analytic formalism firmly entrenched in 



QCD. In fact, some authors 0, [3| have studied Y(4T40) via QCD sum rules soon after its observation, 
however, they arrived at different conclusions basing on the D*D* molecular picture. On the other hand, it 
could not be readily excluded for D S D S or D* S D S as possible molecular configuration for F(4140) without 
explicit dynamics calculations. Catalyzed by the above reasons, we devote to calculate the spectra of 
the (Qs)(*\Qs)(*^ molecular states through QCD sum rules, to see whether F(4140) can be figured as a 
molecular state. In our approach, the masses for D S D S , D*D S , D*D*, B S B S , B*B S , and B*B* molecular 
states are gained. In addition, to improve on the accuracy of QCD sum rule analysis for V(4140), the m 2 s 
order and (g 3 G 3 ) contributions are included in OPE side. 

The paper is organized as follows. In Sec. [Til QCD sum rules for the molecular states are introduced, and 
both the phcnomcnological representation and QCD side are derived, followed by the numerical analysis 
to extract the hadronic masses in Sec. IIIII Section IIVI is a brief summary. 



II. (Qs) w (Qs) w QCD SUM RULES 

The QCD sum rule attempts to link the hadron phenomenology with the interactions of quarks and 
gluons, which contains three main ingredients: an approximate description of the correlator in terms of 
intermediate states through the dispersion relation, a description of the same correlator in terms of QCD 
degrees of freedom via an OPE, and a procedure for matching these two descriptions and extracting the 
parameters that characterize the hadronic state of interest. 
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A. the molecular state QCD sum rule 

In the molecular pictures, following forms of currents can be constructed for (Qs)(*\Qs)(*^ states, with 

•?'««) (<3«) = ( E aij5Qa){Qbij5Sb), 

for one type of hadrons, and 

J (Qsy(Qs) = (^7 AI 0a)(0^75S 6 ), 

for another type, where a and b are color indices. 

For the former case, the starting point is the two-point correlator 

n( q 2 ) = i [ d 4 ^-(o|r[j(x)j+(o)]|o). (l) 



The correlator can be phenomenologically expressed as a dispersion integral over a physical spectral func- 
tion 

tt, 2x 1 [°° , ImnP hcn ( S ) 

n(<? ) = ; r9 — o + - / ds r>~^ + subtractions, (2) 

Mjj-q 2 tt J S0 s-q 2 

where Mh denotes the mass of the hadronic resonance, and Xh gives the coupling of the current to the 
hadron (Q\j\H) = Xh- In the OPE side, the correlator can be written in terms of a dispersion relation as 

roc n 0PE (s) 

n(q 2 ) = / ds p - p., (3) 

J(2m Q +2m s ) 2 s ~ Q 

where the spectral density is given by the imaginary part of the correlator 

p 0PE ( S ) = ilmn 0PE ( S ). (4) 

TT 

After equating the two sides, assuming quark-hadron duality, and making a Borel transform, the sum rule 
can be written as 

X 2 He -Ml/M> = r dsf) OPE {s)e - s /M^ (5) 

J (2m Q +2m s ) 2 

To eliminate the hadronic coupling constant Xh, one reckons the ratio of derivative of the sum rule and 
itself, and then yields 



Ml = I" dsp 0PE se-^/ I" dsp ^e-^\ (6) 

J(2m Q +2m 3 ) 2 J(2m Q +2m s ) 2 

For the latter case, one starts from the two-point correlator 

IF"V) = i J #xe*- x {Q\T\j»{x)j v+ (W>)- (7) 
Lorentz covariancc implies that the two-point correlation function can be generally parameterized as 

ipV) = (^f - 5 nn (1 V) + n(°'(g 2 ). (8) 

q z q z 

The part of the correlator proportional to will be chosen to extract the mass sum rule here. In 
phenomenology, Il^^g 2 ) can be expressed as a dispersion integral over a physical spectral function 



„ m , 2 , [A«] 2 1 f°° J ImllWP 11011 ^) 

n (1) ((j 2 ) = J „ + - / ds — + subtractions, 

M H -q 2 tt J sa s-q 2 



(9) 
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where Mh denotes the mass of the hadronic resonance. In the OPE side, U^{q 2 ) can be written in terms 
of a dispersion relation as 

too „OPE M 

n«(<z 2 )=/ ds P - W (10) 

J(2in Q +2m,) 2 s 1 

where the spectral density is given by 

p 0PE (s) = ilmn( 1 )(s). (11) 

After equating the two sides, assuming quark-hadron duality, and making a Borel transform, the sum rule 
can be written as 

[A (l) ]2e -Mi/M 2 = / d S p 0FE (s)e- s / M2 . (12) 

J (2m Q +2m B ) 2 

To eliminate the hadronic coupling constant one reckons the ratio of derivative of the sum rule and 
itself, and then yields 

M| = P d S p°™ se -s,M> ! f S0 dgp OPE e ~ s/ M>_ (ig) 

J(2mQ+2m a ) 2 i(2m c +2m s ) s 



B. spectral densities 



Calculating the OPE side, one works at leading order in a s and considers condensates up to dimension 
six with the similar techniques in Refs. 15|, |l6|. The s quark is dealt as a light one and the diagrams 
are considered up to order m 2 . To keep the heavy-quark mass finite, one uses the momentum-space 
expression for the heavy-quark propagator. One calculates the light-quark part of the correlation function 
in the coordinate space, which is then Fourier-transformed to the momentum space in D dimension. The 
resulting light-quark part is combined with the heavy-quark part before it is dimensionally regularized 
at D = 4. For the heavy-quark propagator with two and three gluons attached, the momentum-space 
expressions given in Ref. 17| are used. After some tedious calculations, finally with 



p OPE (s) = pPert (s) + p{ ss) {s) + p{ ssf {s) + p{9 s«. G s) {s) + p (M {s) + p {g 3 G 3 ) 
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for (Qs)*(Qs), and 

p 0PE ( S ) = pP ert (s) +p^(-s) +P {Ss)2 (s) +p &°-G*)(s) +p^ G2 \ S )+p^ G ^{s), 
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for (Qs)*(Qs)*. The integration limits are given by a m in = (1 — Jl-4mg/s)/2, a 
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III. NUMERICAL ANALYSIS 



In this part, the sum rules ([H]) and (|13p will be numerically analyzed. The input values are taken as m c — 
1.23 GeV, m b = 4.20 GeV, m s = 0.13 GeV, (qq) = -(0.23) 3 GeV 3 , (ss) = 0.8 (qq), (gsa ■ Gs) = m\ (ss), 
m\ = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 , and (g 3 G 3 ) = 0.045 GeV 6 . Complying with the standard procedure 
of sum rule analysis, the threshold sq and Borel parameter M 2 are varied to find the optimal stability 
window, in which the perturbative contribution should be larger than the condensate contributions while 
the pole contribution larger than continuum contribution. Thus, the regions of Sq and M 2 are taken 
as ^/sE = 4.3 - 4.5 GeV, M 2 = 3.5 - 4.5 GeV 2 for D s D Sl ^ = 4.5 - 4.7 GeV, M 2 = 3.5 ~ 
4.5 GeV 2 for D* S D S , ^ = 4.6 - 4.8 GeV, M 2 = 3.5 - 4.5 GeV 2 for £>*£>*, = 11.1 - 11.3 GeV, 
M 2 = 9.5 - 11.0 GeV 2 for B S B S , ^ = H-l ~ H-3 GeV, M 2 = 9.5 - 11.0 GeV 2 for B* S B S , and 
^ = 11.2 - 11.4 GeV, M 2 = 9.5 - 11.0 GeV 2 for B*B*, respectively. The corresponding Borel 
curves are exhibited in Figs. 1-3. Ultimately, we obtain the mass values: 3.91 ± 0.10 GeV for D S D S , 
4.01 ± 0.10 GeV for D*D S , 4.13 ± 0.10 GeV for D* S D* S , 10.70 ± 0.10 GeV for B S B S , 10.71 ± 0.11 GeV for 
B*B S , and 10.80± 0.10 GeV for B*B*. It is worth noting that uncertainty in our results are merely owing 
to the sum rule windows (variation of the threshold sq and Borel parameter Af 2 ), not involving the ones 
from the variation of quark masses and QCD parameters. 
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FIG. 1: The dependence on M 2 for the masses of D S D S and B B B S from sum rule (|6]). The continuum thresholds 
are taken as y^So = 4.3 ~ 4.5 GeV and y^io = 11.1 ~ 11.3 GeV, respectively. 
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FIG. 2: The dependence on M 2 for the masses of D^D a and BgB s from sum rule (|13p . The continuum thresholds 
are taken as y^so = 4.5 ~ 4.7 GeV and ^/sq = 11.1 ~ 11.3 GeV, respectively. 



4.6GeV 
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FIG. 3: The dependence on M 2 for the masses of D* S D* S and B^Bl from sum rule ([6]). The continuum thresholds 
are taken as %/so = 4.6 ~ 4.8 GeV and t/so = 11.2 ~ 11.4 GeV, respectively. 

IV. SUMMARY 

In summary, the QCD sum rules have been employed to compute the masses of (Qs)W(Qs)W, including 
the contributions of the operators up to dimension six in OPE. For the charmed molecular states, we have 
got M DsBs = 3.91 ±0.10 GeV, M D , Ds = 4.01 ±0.10 GeV, and M D , B , = 4.13 ±0.10 GeV. The numerical 
values for D S D S and D*D S are lower than the mass of V(4140), 4143.0 ± 2.9 ± 1.2 MeV, whereas, the 
one for D* S D* S is well compatible with the experimental data, which supports the D*D* configuration for 
F(4140). Additionally, we have extracted M BsSs = 10.70 ±0.10 GeV, M BjSa = 10.71 ±0.11 GeV, and 
M B ,g, = 10.80 ± 0.10 GeV for the bottom molecular states. Altogether, all these theoretical results are 
looking forward to further experimental identification. 
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